Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences

Buruli ulcer (BU) is a neglected tropical skin disease caused by infection with Mycobacterium ulcerans. The mode of its transmission into the skin has not been fully clarified, but it is believed to involve mechanical transfer, perhaps by insect bites. M. ulcerans secretes a toxin called mycolactone which causes necrosis of the skin and formation of large skin lesions. Understanding the development of these skin lesions over time is crucial to understand the pathogenesis of BU. We formulate a mathematical model to describe the space-time dynamics of M. ulcerans bacteria, mycolactone and skin cell densities as well as fibrin deposition. We investigate the relative contributions to ulcer formation of direct cytotoxicity of mycolactone and ischaemia following the loss of coagulation control. Our results are consistent with existing biological evidence that higher mycolactone concentrations accelerate cell death. Furthermore, the model indicates that the rate of wound expansion is primarily controlled by the diffusion and net birth rate of bacteria, which sets the overall scaling for the speed of wound expansion, along with two dimensionless constants that define the shape of the trajectory of wound edge position against time, namely the ratio of the mycolactone diffusion coefficient to the bacterial diffusion coefficient, and the ratio of the mycolactone degradation rate to the bacterial net birth rate. When mycolactone diffuses significantly faster than bacteria spread or when mycolactone degrades more slowly than bacteria reproduce, the initial development of the wound is faster. In contrast, the long-term wound expansion speed appears to be relatively unchanged across a wide range of the dimensionless parameters.

Epithelial monolayers are the foundation of many mammalian organs and have a significant impact on numerous biological processes, including the secretion of cytokines, absorption of waste products and barrier function. During organ development epithelia are subject to many external forces due to the growth of the surrounding tissues. They exhibit complex responses to such forces, with deformation occurring on short timescales (seconds), and relaxation and remodelling occurring on longer timescales (minutes or hours). Existing mathematical and computational models do not typically account for subcellular remodelling, assuming instead that the response to an applied force or strain acts on a single timescale. In this paper we extend an off-lattice, cell-centre modelling framework to account for subcellular and tissue remodelling by introducing a Dynamic Reference Frame (DRF). The DRF provides a discrete, cell-based analogue of morpho-elasticity: it decouples the evolution of each cell's mechanical reference configuration from its instantaneous deformation, providing a phenomenological description of subcellular remodelling processes - cytoskeletal reorganisation, myosin turnover and adhesion bond remodelling - at the cellular scale. Using this extended multicellular model, we reproduce multiscale responses to external mechanical forces, including creep and stress relaxation experiments, and demonstrate that the history of applied deformation influences tissue recovery upon release. Additionally we identify and quantify where and how these multiscale responses occur. The resulting framework allows for more detailed descriptions and analyses of the development and function of biological tissues.

In many biological, biomedical and industrial systems, particles are transported via fluids through confined networks, in which clogging can disrupt function. However, we lack a predictive theoretical framework that couples particle transport, suspension rheology and network flow resistance. Here, we develop a model and solution algorithm for particle suspension flow in networks based on vessel-level continuum modelling and particle distribution at nodes connecting vessels. We apply the model to study transport of dense particle suspensions in minimal and physiological biological networks. A key feature of our model is the coupling between particle volume fraction and particle flux: in line with the physics of dense suspensions, each network branch, or vessel, possesses a local carrying capacity for particle transport at an intermediate particle fraction between zero and the maximum packing fraction. If this flux capacity is reached, the vessel becomes flux-limited and particles can accumulate in upstream branches, causing them to enter a high-particle-fraction, high-resistance state that we refer to as 'clogged'. We show that these vessel flux limitations lead to network-level redistribution of particles, which can cause widespread clogging and emergent network-scale heterogeneity. By varying network topology, we find that in some regimes increasing network connectivity does not improve transport: paradoxically, additional pathways can promote clogging and reduce network-level particle flux, analogous to classic results in traffic flow networks. Our results provide a minimal mechanistic framework that links suspension physics, network topology, and transport failure in complex flow networks.

The study of tissue dynamics has been stimulated during the last decades thanks to the use of quantitative descriptions, with the development of several theoretical and computational frameworks, many of them revolving around the notion of reaction-diffusion systems, eventually with additional structure variables beyond time and space. The use of structure variables can accommodate phenotypic traits. In this work, we study a family of competition models, where a given population depends on a resource (e.g. oxygen) and several populations are competing for it. Our quantitative description incorporates phenotypic traits and heterogeneity at the level of cell cycle variations, which influence replication rates via oxygen consumption. This enables us to replicate the fitness of specific subpopulations to environmental conditions (e.g. oxygen shortage or external influences). Using numerical simulations, we show that such models display dynamical pattern formation in the form of coupled travelling wave profiles that expand or retreat at the same wave speed. The full theoretical analysis of such dynamics is quite involved; to circumvent this difficulty, we introduce a quasi-stationary approximation for the resource dynamics. We find that this approximation can reproduce the overall behaviour very accurately, with the additional benefit of allowing theoretical treatment of the reduced model. In this way, we provide estimates on the wave speed which are numerically shown to be robust across a wide range of macroscopic parameters of the full model. The wave speeds are thus found to depend strongly on the proliferation rate of the fittest population, resembling a winner-takes-all dynamics.

BackgroundPlatelet activation in flowing blood is a multiscale process in which vessel-scale hemodynamics, red blood cell (RBC) mechanics, adhesive receptor interactions, and intracellular signalling jointly determine thrombotic risk. Individual components are well studied, but a single reduced description that carries each explicitly from vessel-scale flow to mechanosensitive calcium entry, with dimensionally consistent couplings, remains uncommon. ObjectivesWe develop and analyse a reduced, six-module mechanobiological framework for platelet priming spanning the cascade from hemodynamic shear to mechanosensitive calcium entry, and we delineate which elements are supported by existing evidence and which are new, testable hypotheses. MethodsThe framework comprises six coupled modules: (I) hemodynamic forcing from the incompressible Navier-Stokes equations, with an objective principal-strain-rate measure for extensional flow; (II) RBC-mediated platelet margination and near-wall delivery, closed by a near-wall arrival flux; (III) von Willebrand factor (VWF) activation with a bounded kernel and glycoprotein Ib (GPIb) catch-slip capture, resolved through an explicit contact area and a bond-dependent mobility that progressively immobilises wall-interacting platelets; (IV) a single-load membrane-stimulus formulation; (V) mechanosensitive gating and a dimensionally consistent cytosol-store calcium model with extracellular influx; and (VI) a phenomenological mechanical-memory state. We formally derive that the single-platelet stochastic dynamics and the continuum population balance form a Fokker-Planck pair, with the spatially varying diffusivity handled by an explicit drift correction. ResultsThe framework yields a family of mechanochemical dimensionless groups delineating priming regimes. Its central prediction is reformulated as a falsifiable, history-sensitive signature: in a conditioning-test protocol, a low-tension conditioning block charges the memory state, and a fixed sub-threshold test pulse then reports a delay-dependent calcium facilitation that decays on the memory time{tau} m and is distinguishable from no-memory gating, channel adaptation, and residual-calcium priming. We show explicitly that the previously proposed pulsatile-versus-monotone contrast is a nonlinear convexity/thresholding effect of the gating nonlinearity--its difference-in-differences is approximately zero-- and is therefore not a valid test of memory; the conditioning-test signature is. A second prediction links RBC stiffening to reduced near-wall delivery and captured-platelet calcium response, upstream of intrinsic platelet signalling. ConclusionsThe framework provides a dimensionally consistent, mechanistically grounded and hypothesis-generating description linking hemodynamic forcing to mechanosensitive calcium entry. It demonstrates how history-dependent platelet priming may arise from a phenomenological sensitisation state and proposes a conditioning-test protocol for comparison against adhesive, channel and intracellular-store persistence. The framework is calibratable rather than validated, and the quantitative outputs shown use representative uncalibrated parameters.

We develop a general geometrical framework for growth processes in which the growth rate depends exponentially on the number of active growth-supporting components. The framework separates the fraction of active components from the size of the accessible component space, allowing qualitatively different growth laws to arise from the same underlying construction. Bacterial population growth provides an example of a fixed-size component space controlled primarily by activation and deactivation, whereas examples such as human population growth or speciation provide an expanding component space in which interactions between individuals promote further growth. This motivates a Replicative geometrical law, dX/dt= X exp({beta} X^{gamma}), whose exponent characterizes the scaling of the effective growth-supporting component space and is tested against historical human-population regimes and eukaryotic lineage-through-time data. For innovation indicators such as patents and arXiv manuscript submissions, where an existing item need not itself reproduce, the same construction motivates a Direct Geometrical production law without the prefactor X. The framework therefore provides a common description of fixed and expanding component spaces while distinguishing Replicative from Direct Geometrical growth.

Structural network representations of metastatic dissemination typically focus on static topology without resolving transport dynamics, relaxation timescales, or steady-state behaviour. Here, we formulate a discrete Markovian transport model on a directed higher-order network with transition rates derived from qualitative clinical affinity classes. By constructing a non-Hermitian row-stochastic transfer operator, we characterise the relaxation dynamics through its spectral decomposition. The system exhibits a fast-mixing regime characterised by a spectral gap of {gamma} {approx} 0.67, corresponding to a characteristic relaxation timescale of {tau} {approx} 1.49 discrete steps, with the influence of the primary tumour origin progressively attenuated during dissemination. Convergence towards a non-equilibrium steady state (NESS) is accompanied by a reduction in Shannon entropy, concentrating probability mass within specific topological sinks. This spectral relaxation delineates two distinct dynamical regimes: early transient dissemination (n < {tau}), dominated by local organ-specific transition probabilities (organotropism), and the asymptotic regime (n > {tau}), determined increasingly by the global transport architecture of the network. Comparison with independent clinical and autopsy observations across 21 primary tumours and 23 target organs indicates that the predicted stationary distribution is consistent with the observed hierarchy of metastatic organ involvement.

Biological processes, from molecules diffusing to their regulatory destinations to whales pursuing food or mates, are widely exploratory. Such behaviors are effective when some verifiable functional state or outcome can be reached regardless of initial conditions. In these cases, systems repeatedly undergo distinct and abortive trajectories; the process only ends when the system finds 'right' outcomes. This dominance of the terminal condition (and indifference to the initial state) provokes a very different perspective than the conventional 'dynamical systems' framework that has been a centerpiece of the quantitative sciences for centuries. We hypothesize that many of these problems defy the initial-condition driven or gradients on landscapes so useful in problems ranging from mechanics to electrodynamics to chemical kinetics to mass and heat transport. We examine several mathematical frameworks that capture and unify key aspects of exploratory dynamics. One powerful way of thinking of such processes is the geometric distribution, where repeated failures are punctuated by a successful trajectory. We develop intuitions for how random walks with resets accelerate search processes. We highlight fresh and surprising behaviors of search under drift, cues, and checkpoints. Last, appreciating the probability of trajectories conditioned on satisfying macroscopic final outcomes reveals a language for the potency of variation sculpted by selection in their broadest forms. This can give the fictitious appearance of the future making itself known in the present, but we view this as conceptually similar to the way 'fictitious forces' arise in non-inertial reference frames. These approaches stress unity, open questions, and applications across a range of biological phenomena.

Healing of the gastrointestinal (GI) tract following surgical interventions is governed by coupled mechanobiological processes and may lead to severe complications such as perforation or stricture formation. A first comprehensive mechanobiological modeling framework of GI healing is herein proposed. The framework addresses two clinically relevant applications of bioprinting: endoscopic treatment of GI tissue defects and bioprinting-assisted surgical anastomosis. Tissue remodeling is based upon a hyperelastic material model with fiber dispersion coupling the spatiotemporal reaction diffusion dynamics of fibroblasts and cytokines with collagen deposition, reorientation, and plastic deformation. An active stress formulation drives fibroblasts and myofibroblast contraction while a multiplicative decomposition approach is adopted for collagen remodeling as a permanent change. The governing equations are discretised in time with an implicit Euler scheme, finite element in space and solved by a staggered algorithm implemented in the open-source platform FEniCS. Numerical simulations show that the initial collagen fibers content and orientation strongly influence fibroblast infiltration and scar contracture formation. The study quantifies the advantage of having a collagen bioprinted architecture aligned with tissue fibers, to reduce structural remodeling distortion. The model also predicts a strong correlation between collagen deposition and luminal burst pressure during the recovery of mechanical integrity in GI anastomosis.

Large shear deformation in cartilage can lead to tissue failure and injury, which are often precursors to diseases such as post-traumatic osteoarthritis. Under shear, cartilage strain-stiffens prior to failure, which makes stiffening an important indicator of changes to the loading of the extracellular matrix that gives rise to this tissue's properties. While a variety of stiffening mechanisms ranging from uncrimping to fiber-alignment have been explored, whether such mechanisms apply to the depth dependent shear strain stiffening in articular cartilage remains poorly understood. Here, we present experiments and theory that support the idea that the strain-stiffening behavior is caused by a mechanical phase transition from a random isotropic network whose mechanics arise primarily from bending-dominated modes to one where the mechanics arise from stretching-dominated modes. Indeed, consistent with network model predictions, we find that stress-strain data for all depths can be collapsed onto a single universal curve that diverges at a depth-dependent critical strain. These results show that increasing the local concentration of collagen, which constitutes the primary load-bearing network in cartilage, decreases the critical strain at which the strain-stiffening transition occurs, while increasing the low-strain elastic modulus. Furthermore, we demonstrate that bulk compression leads to a higher critical strain and decreased low-strain elastic modulus. Network simulations capture this rich behavior, and indicate that under compression the network collapses, which allows for larger strains before the transition from bending-dominated to stretching-dominated mechanics occurs. Collectively, these findings shed light on the roles that tissue composition and axial strain play in the mechanisms that organize how articular cartilage stiffens during the lead-up to injury.

Devireddy et al. [1] previously showed that an enthalpy-based, macroscale model simulated tissue freezing histories are in close agreement with a coupled model incorporating cellular water transport and intracellular ice formation (IIF). That comparison, however, was performed for a limited set of tissue geometries, cooling conditions, and biophysical parameters. Here, we revisit that conclusion by reimplementing the published model from its governing equations and systematically expanding the parameter space to determine when microscale processes alter the predicted macroscale freezing response. We incorporate both surface catalyzed nucleation (SCN) and volume catalyzed nucleation (VCN) of intracellular ice with cellular osmotic dehydration and examine a broad range of tissue dimensions, cooling rates, convective boundary conditions, membrane permeabilities, cell sizes, activation energies, intracellular water fractions, and nucleation rates. Linearization of the cellular water transport equation about osmotic equilibrium yields an osmotic relaxation timethat can be compared with the local residence time of tissue undergoing phase change. Their ratio defines a dimensionless group, that organizes the transition between regimes in which the uncoupled enthalpy formulation and the coupled micro and macroscale formulation produce similar or substantially different results. Across the one-dimensional parameter space examined here, the two models remain within a prescribed thermal-error bound for the dimensionless group below approximately 0.1, whereas the discrepancy increases systematically above this threshold. Model divergence also depends on intracellular ice nucleation, i.e., when nucleation is sufficiently slow, cellular water transport can become rate limiting for latent heat release. Conversely, sufficiently rapid intracellular ice formation reduces the macroscopic discrepancy by providing an additional pathway for latent-heat release. Multidimensional simulations further show that geometry alone can move a nominally conventional freezing protocol into the coupled regime. In particular, multidirectional cooling can reduce the local residence time sufficiently to produce substantial differences in predicted frozen volume even when the corresponding thermal histories appear similar. Finally, we formulate an inverse approach for estimating IIF parameters from tissue scale thermal histories. Synthetic simulations indicate that parameter recovery is feasible within a finite range of cooling conditions, providing a testable strategy for future experimental measurements. Together, these results provide a mechanistic criterion for deciding when microscale biophysical processes must be retained in macroscale tissue-freezing models.

The relatively long flagella of spermatozoa from insects, birds, and octopuses display double waves, characterized by two superimposed helical waves. The prevalance of these highly organized waveforms across diverse taxa and distinct flagellar architectures hints at shared underlying physics, motivating a model of the flagellum as an elastic filament immersed in a viscous fluid, actively driven by a single set of internal bending moment waves. Simulations of a clamped filament show that it can buckle under its own activity into whirling and flapping states. A multiple-scales analysis of the elastohydrodynamic equations reveals how nonlinear interactions between fast undulations generate an effective compression driving buckling, and connects wave-driven buckling to classical follower-force instabilities. Extending the model to a swimming spermatozoon, the same instability produces double waves. Parameter estimates across species show that most observed double waves lie within the regime where buckling is permitted, supporting self-buckling as a generic physical mechanism for double waves.

We investigate the dynamics of the spread of the infection due to a phage- during the range expansion of an E coli colony. We find a rich panoply of spatio-temporal dynamics at large colony-level length scales, arising due to the interplay of mechanical forces, growth and viral infection dynamics at smaller cellular length scales. Strikingly, the phages "surf" the front of the growing bacterial colony by hitchhiking on E coli that are advected due to the growth of bacteria in the colony, resulting in an anisotropic spread of the phages. We identify microscopic processes - of the phage release during lysis (occurring on short millisecond timescales at length scales close to that of a single bacterium) and local nematic alignment of the rod-like E coli bacteria (occurring on timescales comparable to the bacterial growth rate and the length scales of a few bacterial cells) - that enhance the advective effects driving the hitchhiking behavior. Combining our experiments with a mathematical model, we explain the multiple spatio-temporal dynamical regimes - from the coexistence of uninfected bacteria, resistant (lysogenic) cells and phages, to fixed points where the entire population turns resistant or remains uninfected. Our work broaches a new frontier in the investigation of eco-evolutionary processes arising from the mechanical and ecological forces due to the interactions between viruses and their hosts.

Mechanical properties of biological tissues, driven by passive and active forces, play a vital role in several processes ranging from development to cancer metastasis. However, the dynamical responses of cells in tissues, subject to mechanical deformations such as shear and the associated rheological properties, are not well characterized. Here, we use three-dimensional agent-based models for normal and cancer tissues to investigate their responses to simple shear as a function of cell stiffness and stochastic active forces. In the normal epithelium, with uniform strength of active force, the yield stress as a function of shear rate follows the Herschel-Bulkley form over a range of cell volume fraction. Strikingly, the shear rate dependence and the elasticity-dependent changes in the yield stress fall on master curves upon suitable scaling. To model cancer-like behavior, a certain fraction (Np) of cells was chosen to have enhanced activity and decreased stiffness. As Np increases, the extent of collective cell movement decreases, transitioning from affine (collective) to non-affine (individualistic) movement, a finding that is in accord with imaging experiments. Simulations of a model of a stiff solid tumor, with radius Rs embedded in normal tissue, show that as Rs increases, the yield stress increases. Interestingly, the cells migrate collectively as Rs increases. A Gaussian Mixture Model (GMM) and a mean field theory quantitatively account for the simulation as well as experimental results on cancerous, non-cancerous, and a mixture of these two types. The combined theoretical and experimental study establishes that heterogeneity in stiffness and activity determines non-affine movements in normal and cancer tissues.

We investigate the mechanical principles underlying fruiting body morphogenesis in Dictyostelium discoideum. Quantitative shape analysis based on the Young--Laplace law, together with AFM indentation measurements, indicate surface tension as the dominant tissue-scale force acting on the culminating fruiting body. Based on this observation, we construct a hydrodynamic phase-field model with tunable surface and interfacial tensions, and analyze its behavior numerically. Our results show that, once a stalk begins to form, the elevation of the cell mass arises naturally through a dewetting process. Through quantitative comparisons with experimental measurements, we identify the mechanical conditions required for detachment from the substrate and for establishment of the characteristic morphology of the culminating fruiting body. Together, our model analysis highlights the importance of stalk-tip elongation and tissue-scale surface and interfacial tensions in the construction of large-scale three-dimensional tissues.

Much of the literature and intuition associated with mathematical epidemiology is driven by deterministic models, which are a reasonable assumption when the population size is large. Stochastic models, especially individual based models, are however considered vital when dealing with small population sizes, especially at times of invasion or extinction. The overwhelming majority of these models (both deterministic and stochastic) assume that the underlying parameters are fixed (or follow a regular seasonal pattern). Here, we consider an analytic framework for dealing with randomly varying parameters through the use of stochastic differential equations - thereby capturing the action of external noisy processes such as weather. In particular, we focus on when the transmission rate, {beta}, varies as the solution to a Cox-Ingersoll-Ross Model, such that {beta} is gamma distributed with autocorrelation. We consider the impact of this parameter variation on a stochastic version of the Susceptible-Infected-Recovered model, and for this 'double-stochastic' model show through simulation and analytical results that exogenous noise increases the impact of stochasticity, potentially leading to more early extinctions, wider variations in the number of cases at equilibrium, but that early growth rate can be faster or slower depending on the precise parameters.

A fundamental biological process of transcription occurs in the cell nucleus, which is a complex medium that also contains multiple heterogeneous structures known as biomolecular condensates. Interestingly, some of these condensates contain mRNA molecules in addition to proteins, suggesting an important cellular role in transcription that is not yet well understood. In this work, we develop a minimal theoretical framework for quantitative investigation of the role of reversible mRNA condensation in transcription. Our discrete-state stochastic approach accounts for the most relevant processes, allowing us to explicitly evaluate the properties of the system and clarify the effects of condensation. Analytical calculations supported by computer simulations suggest that reversible mRNA condensation influences the transcription processes by maintaining a constant level of free mRNA in the nucleoplasm while lowering the degree of stochastic noise and increasing the robustness against external perturbations. Physicochemical arguments are presented to explain these observations. The proposed theoretical framework elucidates important microscopic aspects of transcription, providing a convenient quantitative tool for investigating complex biological phenomena.

Purpose: We introduce the Temporal Index of Divergent Evidence (TIDE), a metric that evaluates how the topic center of gravity of a corresponding corpus of literature aligns with, or shifts relative to, a recommendation's semantics over time. Methods: TIDE computes annual mean cosine similarity using MedCPT embedding vectors between a recommendation and corresponding PubMed records binned by year (2000-2024), testing for monotonic temporal trend (Kendall {tau}). Three recommendations with expected divergent, stable, and convergent trajectories were examined: two-stage exchange for periprosthetic joint infection (PJI), emergent fasciotomy for acute compartment syndrome (ACS), and non-operative management of distal radius fractures (DRF) in adults [&ge;] 65 years. The signal was characterized by an additive per-token attribution, an isolated vs. in-context analysis, and an adversarial minimal-pair test. Results: PJI diverged ({tau} = -0.84, year-shuffle p = 0.001), DRF converged ({tau} = +0.42, p = 0.005), and ACS showed no detectable trend ({tau} = -0.07, p = 0.66). The two significant directions (PJI, DRF) held under matched input [CLS] pooling, while ACS remained near zero. Publication volume baselines did not distinguish the three trajectories (publication count was universally positive, {tau} = +0.77 to +0.97). A reversed-meaning minimal-pair control showed little cosine sensitivity to reversal. Paper-level TF-IDF recovered the qualitative three-case pattern. Conclusion: Embedding-based drift analysis is a proof-of-concept candidate pre-triage signal that tracks topical alignment, not evidentiary direction; superiority over TF-IDF remains unproven. The present study is hypothesis-generating, showing potential use of embedding representations as an analytical unit rather than simply for retrieval.

Tumor spheroid mechanics arise from both the mechanical properties of individual cells and the adhesive interactions that organize them into tissues. The relative contribution of these two factors to the bulk mechanical behavior, however, remains difficult to disentangle experimentally. Here, we develop a computational model of micropipette aspiration to compare the mechanical response of isolated cells and multicellular spheroids within a common computational framework. By independently varying single-cell stiffness and cell-cell adhesion, we quantify their effects on aspiration dynamics, effective elastic modulus, and viscoelastic relaxation. Our results show that increasing single-cell stiffness substantially alters the mechanics of isolated cells but has limited influence on the effective elastic modulus of multicellular spheroids. In contrast, changes in cell-cell adhesion produce pronounced effects on spheroid effective elastic modulus. Nevertheless, both parameters increase the retardation time governing the transition from the initial elastic response to long-time viscous deformation. These findings suggest that multicellular elasticity is governed primarily by intercellular mechanical coupling, whereas the dynamical response to applied stress depends jointly on cell-scale mechanics and cell-cell adhesion.

Contact inhibition of proliferation (CIP) enables epithelial tissues to self-regulate growth and maintain tissue homeostasis. However, how cell-level mechanical contact, tissue-scale structural order, and proliferation kinetics interplay remains a fundamental open question in living matter physics. Here, we present a particle-based model of a confluent epithelial monolayer governed by overdamped dynamics, where individual cells interact via a two-dimensional hard core- soft shoulder potential. By comparing structural evolution during quasistatic densification with previously reported experimental division kinetics, we find that the dynamics of proliferation arrest mimics the onset of direct steric contacts between the hard cores of the shell. Identifying hard core contacts as the physical driver of CIP, we couple our mechanical model with a stochastic Monte Carlo division scheme in which the instantaneous division rate decreases to zero from an intrinsic value as the number of hard core contact increases to six from zero. We demonstrate that for high intrinsic division rates, the cellular densification outpaces mechanical relaxation. This kinetic mismatch drives premature hard-core contact formation, shifts the onset of jamming and contact inhibition to lower packing fractions, and induces increasingly disordered transient configurations before the tissue universally converges to a hexagonal close-packed limit. Our model's predicted division kinetics and structural order evolution are consistent with epithelial monolayer experiments, both reported and our own. This minimal physical framework links single-cell steric contact mechanics directly to tissue-scale growth regulation and structural evolution.